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We investigate the entanglement of the ferromagnetic XY model in a random magnetic 
field at zero temperature and in the uniform magnetic field at finite temperatures. We use the 
concurrence to quantify the entanglement. We find that, in the ferromagnetic region of the 
uniform magnetic field h, all the concurrences are generated by the random magnetic field and 
by the thermal fluctuation. In one particular region of h, the next-nearest neighbor concurrence 
is generated by the random field but not at finite temperatures. We also find that the qualitative 
behavior of the maximum point of the entanglement in the random magnetic field depends on 
whether the variance of its distribution function is finite or not. 

KEYWORDS: entanglement, concurrence, entanglement of formation, thermal entanglement, decoherence, 
random field, XY model 



1. Introduction 

We study in the present paper the pairwise entan- 
glement of the ferromagnetic spin-1/2 XY chain in a 
random magnetic field at zero temperature and in the 
uniform magnetic field at finite temperatures. The en- 
tanglement 1 ' 2 is one of the most interesting features of 
quantum mechanics. It has the property of non-locality 
originating in the principle of superposition. One typical 
example of the state which has the property of the en- 
tanglement is the singlet \tp) = l/v%(|01) - |10)). The 
measurement on one particle of this state affects the 
other particle immediately, even if two particles are far 
away from each other. This non-locality of the entangled 
state puzzled many people including Einstein in the early 
times when the quantum mechanics was born; Einstein 
thought that quantum mechanics was an incomplete the- 
ory because of the non-locality. However, violation of the 
Bell inequality, 3 or more generally the CHSH inequal- 
ity, 4 showed that the non-locality is a reality. Quantum 
mechanics is now widely accepted, having explained a 
lot of phenomena which were not explained by classical 
mechanics. 

The quantum information processing 5 such as quan- 
tum teleportation 6 and super dense coding 7 have been 
heavily studied recently. The entanglement is vital for 
implementation of such techniques. In reality, however, 
the entanglement may be destroyed by some decoherence 
effects 5 such as the thermal fluctuation and a random 
magnetic field. Thus, it is important to know how the 
entanglement is affected by thermal and impurity dis- 
turbances. Effects of impurities on the entanglement are 
also interesting from the viewpoint of the relation be- 
tween quantum coherence and impurities. 

In the present paper, we calculate the entanglement 
between two spins of the ferromagnetic isotropic XY 
chain in a random magnetic field as well as at finite tem- 
peratures. To our knowledge, this is the first to study sys- 
tematically the dependence of the entanglement on the 
randomness in a spin system. Li et al. s studied the de- 
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pendence in the one-electron Anderson model in one di- 
mension. The conclusion that the entanglement increases 
due to the randomness in some parameter regions (see 
below) is common to both studies. (After submitting the 
present paper, we noticed a study on the entanglement 
in random quantum spin-5 chains. 9 It, however, is not 
quite related to the present issue; the study concerns a 
scaling law of the entanglement entropy in the random 
singlet phase.) 

There have been a couple of works which studied the 
effects of the temperature on the entanglement of spin 
systems. Arnesen et al. 10 mentioned that the nearest- 
neighbor entanglement of the anti-ferromagnetic Heisen- 
berg chain can be increased by introducing the temper- 
ature in a uniform magnetic field. Similar work has been 
done by Nielsen 11 on the two-spin Hcisenberg model. Os- 
borne and Nielsen 12 studied the nearest-neighbor and the 
next-nearest-neighbor entanglement of the anisotropic 
XY chain and the ferromagnetic transverse Ising chain. 
Although their main interest is in the entanglement 
near the quantum ground-state phase transition, they 
also mentioned calculation of the entanglement at finite 
temperatures. Yano and Nishimori 13, 14 also mentioned 
a finite-temperature calculation of the nearest-neighbor 
entanglement on the anti-ferromagnetic anisotropic XY 
model. Their results for the nearest-neighbor entangle- 
ment are almost the same as ours. Our conclusion that 
the entanglement also increases due to the thermal fluc- 
tuation in some parameter regions (see below) is common 
to the above-mentioned studies. None, however, com- 
pared the entanglement in a random magnetic field and 
that at finite temperatures quantitatively in the same 
model. 

We compute the entanglement between two spins from 
the nearest-neighbor pair to the fifth-neighbor pair. We 
use the concurrence 15 to quantify the pairwise entangle- 
ment. We find the following: 

(1) In general, the entanglement is decreased as the ran- 
domness is increased; 

(2) In the region of the uniform magnetic field h > 1, the 
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entanglement is increased by the random magnetic 
field and by the temperature. In yet another region 
ft. < 1/2, it is increased by the random field but not 
by the temperature. 

(3) Qualitative behavior of the maximum point of the 
entanglement depends on the random magnetic 
field, in particular, whether the variance of the dis- 
tribution function is finite or not. 

(4) The entanglement between two spins at finite tem- 
peratures becomes weaker than that in the random 
magnetic field at zero temperature as the distance 
between the two spins gets greater. 

The present paper is organized as follows. In §2, we intro- 
duce the model and review the outline of computation of 
its correlation functions, which are necessary to quantify 
the concurrence. In §3, we calculate the concurrence for: 
in §3.1, the XY spin chain in the uniform magnetic field 
at zero temperature; in §3.2, the XY spin chain in a ran- 
dom magnetic field at zero temperature; in §3.3, the XY 
spin chain in the uniform magnetic field at finite tem- 
peratures. Finally, we give a summary and discussions in 
§4- 

2. The model and the entanglement 

2. 1 Diagonalization of the model Hamiltonian 

The Hamiltonian of the XY spin chain in a random 
magnetic field is given in the form 

N N 

H = - 4 E + - 2 E ( h + w °i> 

j=l 3 = 1 

(1) 

where J (> 0) is the coupling constant, N is the number 
of the spins, a a (a = x,y,z) are the Pauli matrices, h 
is the uniform magnetic field and {hj} are the random 
magnetic field. We impose the periodic boundary condi- 
tions: 

°n+i =<ri, {a = x,y,z). (2) 

Hereafter, the coupling constant J is set to one. The ran- 
dom magnetic field hj at each site obeys the distribution 
function 

P q , a (hj)~ [a 2 - (1 - q)^ 2 }^ , (3) 

where the parameter q determines the type of the distri- 
bution function and a determines the width of the dis- 
tribution function. In particular, eq. (3) is reduced to a 
Gaussian distribution function as q — > 1. In this case, the 
scale parameter a is its standard deviation. Equation (3) 
is also reduced to a Lorentzian distribution function for 
q = 2. In this case, the scale parameter a is its half width 
at half maximum. The variance of the distribution func- 
tion (3) diverges for q > 5/3 and is finite for q < 5/3. In 
the case of the Lorentzian distribution q = 2, Nishimori 16 
analytically calculated the average one-point correlation 
function and obtained lower bounds of the average two- 
point correlation functions. The results in the paper, 16 
however, are not used in the present paper, since we take 
the random average of the concurrence, which is a non- 
linear function of the one-point and two-point correlation 
functions. 



Author Name 

We diagonalize the Hamiltonian (1) as follows. The 
Hamiltonian can be expressed by the Fermi operators 
a) and a after the Jordan- Wigner transformation. The 
Hamiltonian is then reduced to the quadratic form 

N 

with 

A = 

( -h-fo -\ ••• ±\ \ 
-\ -h-h 2 -§ '■■ 

o -i •• •• ; 

: • . . . . ' • _I 

\ ±! ■■■ -2 -h-h N ) 

(5) 

where we dropped a constant term in the Hamilto- 
nian (4). The signs of the (l,N) and (N,l) elements 
in eq. (5) are negative when the number of the Fermions 
in the system is even and positive when odd. The Her- 
mitian matrix A is diagonalized by a unitary matrix V. 
We thus have 

N 

H = J2^Uu (6) 

i=l 

where the operators c\ and a are given by 

N N 

i=i i=i 

and satisfy the anti-commutation relations {cj, cj} = 5ij 
and {cj, Cj} = 0. 

2. 2 Correlation functions and the two-site density ma- 
trix 

In order to quantify the entanglement, we use the con- 
currence 15 related to the entanglement of formation. 17 
The concurrence between the spins at sites i and j is 
calculated from the two-site density matrix pij as 

Cij — max {0, Ai — A 2 — A 3 — A4} , (8) 

where {A^} 4 =1 are the square roots of the eigenvalues of 
the matrix R = PijPij in non-ascending order, Ai > A 2 > 
A 3 > A 4 with pij = (a? O crj) p* (<j\ ® crj) . The complex 
conjugation is taken in the a z basis. The two-site density 
matrix p^ is defined by 

(Hi = Tr j} .p, (9) 

where Trv- denotes the trace over the degrees of freedom 
except for the sites i and j, and p is the density matrix 
of the whole system: p = e~P H jZ. 

The two-site density matrix can be expanded in terms 
of the identity matrix and the Pauli matrices as 

1 3 

Pij = 1 E Pa P a ? a j ' ( 10 ) 

a, 0=0 
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where of denotes the identity operator on the site i, aj — 
of , of = of and of = of. The coefficients p a p are real 
numbers determined by 



(11) 



Hence, sixteen coefficients are needed to determine the 
two-site density matrix in general. Thanks to the sym- 
metry of the Hamiltonian (1), the number of the coeffi- 
cient is reduced to four; we need (of), (of) and (of cr") 
(a = 1, 3) only. The others are zero. 12 Hence, the two-site 
density matrices of the model take the form 



Pij 



\(lij + (at)cr?®I j + ((T*)I i <i 



(12) 



a=l 



where a 1 = o^, ct 2 = o^, ct 3 = o- z and (of oj) = (ofa 2 ). 
The correlation functions are obtained as follows: 18-20 



G 4 



1 -X x\ 



Gi.i 

Gi 



Gi 
Gi 



3,i 

Gis+\ 
Gi+is+i 



G 



j-la+l 



Gi t i+2 
Gi+lj+2 

Gi 



Gi,j 
Gi+ij 



G 



3-1,3 



(13) 
(14) 

(15) 



*j-l,i+2 

where Gij are given as follows in the three cases: 

i) The XY spin chain in the uniform magnetic field 
at zero temperature. In the thermodynamic limit, 
{Gij} are given by 18-20 

J 1 for h > J, 

M ~ \ -1 + f arccos(-^) for h < J, 



for h > J, 
£7^™[(*-J')a«5COs(-$)] 

for /i < J. 



(16) 



(17) 



ii) The XY spin chain in the uniform magnetic field 
at finite temperatures. In the thermodynamic limit, 
{Gij} are given by 18-20 



Gij 



-Sij + 



2 r 

- \ d^- 

JO 1 



cos(i - j) 



+ exp(— (3( J cos (fi + h)) 
(18) 



iii) The XY spin chain in a random magnetic field at 
zero temperature. In this case, {Gij} are given by 



N G 

Gij=2 , £ i V a V il -5 ij , 



(19) 



where V is the matrix diagonalizing the matrix A 
in eq. (5) and Nq is the number of the Fermions. 
In the ground state, the Fermions are filled in the 
levels with ej < in eq. (6). 



2.3 The concurrence 

Now that the coefficients in eq. (10) have been ob- 
tained, we can evaluate the concurrence. We here define 
the average concurrence as the random and spatial aver- 
age: 



C(r) 



1 N 

i=l 



[Ci. 



i+r\ 



(20) 



where [• • • ] av denotes the random average, Cij denotes 
the concurrence between the sites i and j, and N is the 
number of the sites. In the absence of the random mag- 
netic field, the Hamiltonian possesses the translational 
invariance and the averaging is not necessary. We note 
that the entanglement of formation 17 after the random 
average and the spatial average is always greater than 
that obtained by substitution of the average concurrence 
into the relation between the entanglement of formation 
E and the concurrence C 



E{C)=- 



i + Vi-c 2 



log 2 



i + Vi - c* 2 



Iz^E*^ Iz^E* , (21) 



since E(C) is a concave function of the concurrence C. 
Hereafter, we simply refer to the average concurrence as 
the concurrence. 

In the case of no random magnetic field, we calculated 
the concurrence rigorously in the thermodynamic limit. 
In the case with a random magnetic field, we numerically 
evaluated the sample average of the concurrence (20). For 
all the results below in the random case, the number of 
the sites iV is 500 and the number of the samples is 10000. 
In Fig. 1(a), the next-nearest-neighbor concurrence C(2) 
is plotted for q = 2 (the Lorentzian distribution) with 
error bars at h — 0, 0.5, 1, 1.5, 2, 2.5 and 3, but the er- 
rors are almost invisible. In Fig. 1(b), the next-nearest- 
neighbor concurrence C(2) is plotted for q — 2 and for 
the system size N = 10, 100, 250 and 500 with the ran- 
dom average over 10000 samples. The finite-size effect 
is invisible for N > 100. We hence conclude that 10000 
samples and the system size N = 500 arc substantial. We 
calculated the average concurrence C(r) for 1 < r < 5 
but not all the results are plotted below. 

3. Numerical results 

3. 1 XY spin chain in a uniform magnetic field at zero 
temperature 

We first study the concurrence of the XY spin chain in 
the uniform magnetic field. Figure 2 shows the nearest- 
neighbor concurrence C(l), the next-nearest-neighbor 
concurrence C(2), the third-neighbor concurrence C*(3) 
and the fourth- neighbor concurrence C(4). All the con- 
currences rapidly decrease near h = 1, where the quan- 
tum phase transition occurs, and vanish in the region 
h > 1. In the region h > 1, the ground state is given by 
the tensor product of the one-spin state | fi) as 



|GS> = |Ti>|T2>|Ts>- 



(22) 
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Fig. 1. (Color online) (a) The next- nearest-neighbor concurrence 
C(2) for q = 2 (the Lorentzian distribution) with the scale pa- 
rameter a = 0.3, 0.7, 1. All the lines are plotted as a function of 
the uniform magnetic field h. The random average is taken over 
10000 samples with the system size 500. (b) The next-nearest- 
neighbor concurrence C(2) for q = 2. The lines are plotted with 
the scale parameter a = 0.3 with 10000 samples and for the 
system size N = 10, 100, 250, 500. 



where the state | fj) denotes the eigenstate of the ma- 
trix of satisfying the eigenequation of | Ti) = I Ti)- Since 
there is no superposition involved, the entanglement van- 
ishes in ft, > 1. 

3.2 XY spin chain in a random magnetic field at zero 
temperature 

Next, we study the concurrence in a random magnetic 
field (in addition to the uniform magnetic field h) at 
zero temperature. The random magnetic field obeys the 
distribution function (3); we investigate the cases for q — 
1, 1.35, 5/3, 1.85 and 2. 

The nearest-neighbor concurrence C(l) in all cases be- 
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Fig. 2. (Color online) The concurrence of the XY spin chain in 
a uniform magnetic field as a function of the uniform magnetic 
field h. 
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Fig. 3. (Color online) The nearest-neighbor concurrence C(l) of 
the XY spin chain in a random magnetic field at zero tempera- 
ture; (a) for q = 1 (the Gaussian distribution); (b) for q = 2 (the 
Lorentzian distribution). All the data are plotted as functions of 
the uniform magnetic field. 
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Fig. 4. (Color online) The next-nearest-neighbor concurrence C(2) in a random magnetic field at zero temperature; (a) for q = 1 (the 
Gaussian distribution); (b) for q = 1.35; (c) for q = 5/3; (d) for q = 1.85; (e) for q = 2 (the Lorentzian distribution). All the data are 
plotted as functions of the uniform magnetic field. 



haves similarly. In the region h < 1, the nearest- neighbor 
concurrence in a random magnetic field for each q de- 
creases as the distribution width a is increased. We here 
show in Fig. 3 only the cases q = 1 and q = 2. The reduc- 
tion of the nearest-neighbor concurrence C(l) is greater 
as the scale parameter a is increased. For h > 1, the 
nearest-neighbor concurrence for all q is increased by the 
random magnetic field. That is, the random magnetic 
field increases the quantum correlation. The reason why 



the nearest-neighbor concurrence is increased for h > 1 
may be as follows; the random magnetic field flips some 
of the aligned spins of the ferromagnetic (or classical) 
ground state (22) and thereby the flipped spins and their 
neighboring spins restore the quantum interaction. 

As shown in Fig. 4, the next-nearest neighbor concur- 
rence C(2) for each q is decreased for 0.5 < h < 1 as the 
scale parameter a is increased. On the other hand, it is 
increased for h < 1/2 and ft, > 1 as the randomness a 
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Fig. 5. (Color online) Smaller arrows (blue) indicate spins. Bigger 
arrows (grey) indicate the random magnetic field, (a) The ran- 
dom magnetic field happens to be almost constant over a spatial 
region, (b) The random magnetic field is extremely strong at one 
particular site. 



is increased. This is in contrast to the finite-temperature 
case in the next subsection, where we show that the next- 
nearest-neighbor concurrence C(2) for h < 1/2 is not 
increased at finite temperatures. 

The reason of the increase for h > 1 may be the same 
as the above-mentioned reason for C(l). We have not 
been able to determine decidedly the real reason why 
the next-nearest- neighbor concurrence C(2) is increased 
for h < 1/2. We, however, can consider some situations 
where the next-nearest-neighbor concurrence is increased 
for h < 1/2. Let us first consider the case where the vari- 
ance of the distribution function is finite, i.e. q < 5/3. 
There may arise a situation where the random magnetic 
field is almost constant and equal to the distribution 
width a over a region of considerable length; see Fig. 5(a). 
The uniformness then may generate the concurrence in 
the region. We can give an argument for this specula- 
tion. In Fig. 4(a) and (b), the concurrence without the 
uniform field, h = 0, is generated only for a > 0.5. This 
is consistent with the fact that the concurrence without 
the randomness is zero for h < 0.5; if our speculation is 
correct, the concurrence is generated by an "almost uni- 
form" random field only when the field is greater than 
0.5. 

In the case where the variance of the distribution func- 
tion is infinite, i.e. q > 5/3, we could think of a more 
plausible situation. In this case, a singularly strong ran- 
dom field can appear at a site as illustrated in Fig. 5(b). 
We then can take the Zeeman energy of the site as the 
non-perturbation term and calculate the second-order 
perturbation of the exchange energy. We may end up 
with an effective interaction between the two spins be- 
side the strong magnetic field. In this situation, the next- 
nearest-neighbor concurrence may be restored around 
the strong field. 

We find that the qualitative behavior of the next- 
nearest neighbor concurrence C(2) is different depend- 
ing on whether the variance of the distribution function 
is finite or not. The maximum point of the next-nearest- 
neighbor concurrence for q < 5/3, where the variance of 
the distribution function is finite, shifts to the right as 
the randomness a is increased as shown in Fig. 4(a) and 
(b). In contrast, the maximum point of the next-nearest- 
neighbor concurrence for q > 5/3 in Fig. 4(c)-(e) shifts 
to the left as the randomness a is increased. 
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Fig. 6. (Color online) The third neighbor concurrence C(3) in a 
random magnetic field at zero temperature; (a) for q = 1 (the 
Gaussian distribution); (e) for q = 2 (the Lorentzian distribu- 
tion). All the data are plotted as functions of the uniform mag- 
netic field. 



The third-neighbor concurrence C(3) and the rest, 
C(4) and C(5), behave similarly to the next-nearest- 
neighbor concurrence C(2), only smaller than the next- 
nearest-neighbor concurrence. The maximum point of 
the third-neighbor concurrence and the rest for the cases 
<7 < 5/3 first shift to the left and turn to the right as 
shown in Fig. 6. 

3.3 XY spin chain in a uniform magnetic field at finite 
temperatures 

Third, we investigate the concurrence of the XY spin 
chain in the uniform magnetic field at finite tempera- 
tures; see Fig. 7. The nearest-neighbor concurrence C(l) 
is decreased for h < 1, whereas it is increased for h > 1; 
see Fig. 7(a). The reason why the concurrence is in- 
creased for h > 1 may be that the temperature excites 
entangled states above the ferromagnetic ground state 
with some probability due to thermal fluctuation. Thus, 
the concurrence can have a non-zero value for h > 1 . We 
can hardly see the essential difference between the effects 
of the random magnetic field and the temperature on the 
nearest-neighbor concurrence C(l); compare Fig. 3 and 
Fig. 7(a). 
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Fig. 7. (Color online) The concurrence in a uniform field at fi- 
nite temperatures; (a) The nearest-neighbor concurrence C(l); 
(b) The next-nearest- neighbor concurrence C(2); (c) The third 
neighbor concurrence C(3). 



The next-nearest neighbor concurrence C(2) is de- 
creased for 0.5 < ft. < 1. The increase of the concurrence 
appears only for ft > 1; the concurrence for ft < 0.5 does 
not appear. This is in contrast to the case of the random 
magnetic field shown in the previous subsection, where 
the next-nearest-neighbor concurrence C(2) for ft < 1/2 
is increased by the random field; compare Fig. 4 and 
Fig. 7(b). The difference between the disturbance of the 
random magnetic field and the temperature appears in 
this point. The third-neighbor concurrence C(3) and the 



rest, C(4) and C(5), behave similarly to the next-nearest- 
neighbor concurrence C(2) except for quantitative differ- 
ence. 

Finally, we study the maximum concurrence as a func- 
tion of the scale parameter a or the temperature kT. 
The reduction of the maximum concurrence as a func- 
tion of the scale parameter or the temperature is plot- 
ted in Fig. 8. The nearest-neighbor concurrence C(l) 
decreases for the random magnetic field more rapidly 
than that at finite temperatures in the plotted ranges. In 
the same plot ranges, however, the third-neighbor con- 
currence C(3) in the random magnetic field remains fi- 
nite, whereas the third- neighbor concurrence C(3) at fi- 
nite temperatures almost vanishes for kT > 0.07. We 
thus observe that, as the distance between the two spins 
increases, the concurrence becomes considerably weak 
against the thermal fluctuation. 

4. Summary and discussions 

We have studied the entanglement of the XY spin 
chain in a random magnetic field at zero temperature 
and in a uniform field at finite temperatures. We found 
that: (i) In general, the entanglement is decreased by 
the random magnetic field and the temperature; (ii) The 
entanglement is increased by the random magnetic field 
and the temperature in some parameter regions. That 
is, quantum correlation can be both increased and de- 
creased by the disturbances. In particular, we find that 
the next-nearest- neighbor concurrence C(2) for ft < 1/2 
is increased by the random magnetic field as well as for 
ft > 1 . The increase of the concurrence for ft < 1 /2 does 
not occur for the thermal fluctuation. This is a notable 
difference between the random magnetic field and the 
thermal fluctuation, (iii) The qualitative behavior of the 
maximum point of the concurrence depends on whether 
the variance of the distribution function is finite or not. 
In particular, the maximum point of the concurrence 
shifts to the right when the variance of the distribution 
function is finite, whereas it shifts to the left when the 
variance of the distribution function is infinite, (iv) The 
entanglement between two spins at finite temperatures 
is weaker than that in the random magnetic field at zero 
temperature as the length of two spin is larger. 

In view of implementation of quantum information 
processing, the conclusion (ii) indicates that impurity 
and thermal disturbances are not always destructive to 
the entanglement resources. The conclusion (iii) also in- 
dicates that, under some randomness, we may obtain the 
best quality of the entanglement at a point different from 
the pure case. 

We need further studies to find definite reasons why 
the concurrence is increased by the random magnetic 
field and the thermal fluctuation. It may be also inter- 
esting to compute perturbationally the dependence of 
the concurrence on the distribution width of the random 
magnetic field to see why the behavior of the maximum 
concurrence is different for finite variance and infinite 
variance. 
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Fig. 8. (Color online) (a) The maximum point of the concurrence 
in a uniform magnetic field at finite temperatures as a function of 
the temperature kT; (b) The maximum point of the concurrence 
in a random magnetic field at zero temperature for q = 1 as a 
function of the scale parameter a; (c) The maximum point of the 
concurrence in a random magnetic field at zero temperature for 
q = 2 as a function of the scale parameter a. 



